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Abstract 


In this work the connection between curvature, the emergence of the universe, and 
the process of morphogenesis is investigated in detail. We have come a long way from 
general relativity using differential geometry, which describes the curvature of space- 
time as an emergent phenomena, to the more recent theory of loop quantum gravity 
using non-geometric objects such as holonomy, fluxes, spin, and knotting. However, 
these theories are still only a model for space-time curvature, i.e., they explain the 
coarse-grained description of the structure of space-time and how this can account for 
the scale-variant behaviour and its dynamics. In this work, I will use the term mor- 
phogenesis to describe the emergence of structure in space-time fabric. The different 


levels of scale of phenomena can be described by their fractal analogues. 
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1 Introdcution 


A unified description of gravity is a longstanding problem in physics. Since the early history 
of physics it has been proposed that space-time was discrete, a question that was taken up 
again in the 1980s when physicists started to take the nature of quantum gravity seriously. 
So far, all attempts to find a theory of everything have been plagued by the problem of 
singularities. While these theories are very successful at explaining physics on the largest 
scales of astrophysics, they cannot describe the physics on the smallest scale of atoms and 
elementary particles. One of the latest efforts to find a theory of everything uses a geometric 
approach to finding the quantum mechanics of gravity based on a discrete version of general 
relativity called loop quantum gravity. In these theories, gravity and quantum mechanics are 
unified using non-geometric objects such as holonomy, fluxes, spin, and knotting (see, e.g., 
[1]. However, all of these theories describe the emergent nature of the geometry of space, not 
space as a substance. To explain the full behavior of space, we need to provide an explanation 
of the internal behavior of space. Here I will offer a purely geometric description of this 
emergent process known as morphogenesis that describes the process of space-time developing 
its structure. This article explains the relationship between spacetime curvature, topology 
change, and the fractal nature of space. The key concept of this article is the geometric 
morphogenesis of space-time. This approach leads to a unified understanding of general 
relativity, and the standard model of particle physics as facets of the same phenomena. These 
new ideas replace the current idea of a space-time manifold of differential equations with a 


nonlinear space-time model that is not limited by the smooth and continuous approximation 


of general relativity, but can also be applied to quantum theories of gravity. Unlike many 
other approaches to quantum gravity, which are still using differential geometry in their 
theory, this new geometric approach starts with the assumption of an underlying fractal. 
With a fractal, you have different scales of smoothness, i.e., from coarse to fine. Based on 
the different scales of smoothness, you get a set of differential equations that will allow you 
to describe the degrees of curvature at varying scales. This description is carried out using 


a morphogenetic potential, which gives the geometric curvature of spacetime. 


2 Geometric Morphogensis 
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Figure 1: Morphogenesis of Space-time 


Morphogenesis: It’s defined as the formation of structure in spacetime as a function of 
the geometric morphogenetic potential f, and the external forcing or boundary conditions 
b. The geometric morphogenetic potential is an abstract vector field describing the struc- 
ture formation that emerges. For example, if it is an attractor field, it has a stable space 
of equilibria. If it is a pseudo force field, it has no stable point of equilibria, but instead 
it has a fractal structure with fixed points and saddles. In the former case, the dynamics 


has an eigen spectrum with discrete eigenvalues. In the latter, the eigenvalues are dense on 


the real axis, which is consistent with self-similar fractals. In order to define the geomet- 
ric morphogenetic potential as an abstract vector field, we first identify the n-dimensional 
morphogenetic potential space as 


Qe n 


n times 


We then projectively identify each of the A,,’s along the temporal coordinate, so that 


An becomes identified with R”. The end result is that we have a single infinite dimensional 


vector space whose basis is B,,,; where 7 indexes the dimension of the physical space, e.g.: 


Bni = {point Of(0); t E R} U (2) 
{line l} ,; l+ is a 1-cycle with period Tni} U (3) 
{surface ca ay is a 2-cycle with period o,,;}U (4) 
{volume vi ,; vhs is a 3-cycle with period d,,;} (5) 


The key idea in morphogenesis is that the field acting on a domain gives you not just a 
volume of that domain, but a certain type of shape of that volume. In addition, the field 
also specifies how fast that shape is changing. The general form of morphogenetic fields we 


consider is as follows: 


$ = Ekm k Ak Pkl fkm dV (7) 


tosti 
where fxm, is the geometric morphogenetic potential, p;, is a measure density that we often 
choose to be a Dirac delta function, and Aj", is a measure of the characteristic length-scale 


5 


associated with the geometric morphogenetic potential fkm,- In order to get a non-linear 
equation describing the structure of spacetime, we need to assume that the morphogenetic 


fields are constrained by the incompressibility condition: 


| o dV =1 (8) 


toti 


k 


which means that for a given domain Dý ,,, at every instant in time, the volume associated 


with that domain is a constant, i.e., 


Dina, =U Dea (9) 
l 


Din =U Din (10) 
l 
Volume of domain = | dv=% f dV 
k > 
= X c* = constant (11) 
l 
The temporal scale associated with the domain Df ,, is given by: 
IDE a | 


Ok,ti,t2 = 4| FRE (12) 
a IDE al 


3 Fractals as mathematical tools for Morphogenesis 


Sitot = Ok,ti,t2 Xok a Pitot (13) 


Where m € Ap, k is a topological index which indexes the different spatial domains 
that the geometric morphogenetic field is defined over, and POs. is a component of the 


geometric morphogenetic potential. This model treats the morphogenetic potential as having 


two components, the temporal evolution of the potential, and its spatial distribution. We 
can now extend this definition to n dimensions by taking the composition of n such geometric 
morphogenetic fields. We take the spatial and temporal dimensions of the morphogenetic 
field as one dimension of the metric and a different dimension of the time-space. If m € A, 
is the index for space, and j € A,41 for time-space, the j’th component of f at (x,t) can 
be expressed as follows: 
FE (at) = ous (FG, cb, oh.) 

To allow for inhomogeneities in the morphogenesis process, we also extend f to a function 

that takes in a continuous parameter s from the parameter set S, which then acts as a 


gradient field: 


Sit : (a, O(k,t)) — Jkala, Okt) (14) 


Examples: 


3.1 Inhomogeneity in Geometry 


This is equivalent to taking M as the boundary of spacetime. The external forcings are then 


the initial and final states of the spacetime boundaries: 


F” (s) = (sala), o7"(@) 


With the following boundary conditions 


! 


(0) = Po, ze (1) $ Po, where t < t 


This can be extended to the case of a collection of geometric morphogenetic fields with a 


total of #directions + 1 as follows: 


fig?" = (Eila), hi) 


with the following boundary conditions 


x'(0) = pig, zi (1) = pin where t' < t 


3.2 Self-Organized Criticality in Geometry 


This is equivalent to taking M as the state space of a discrete dynamical system whose 
phase-space is a subset of Euclidean space with non-zero measure EF and a boundary of 


spacetime. The external forcings are then a single point state in the boundary of spacetime: 
m ; ; E 
f™(s) = (zsm (a), pisla), a, a single point £i = (£e; ---, Eeay L— (1 — a)a), where a = g 


The boundary conditions are then trivial: 


f” (0) = (zs) (a), Piola), a, Piz(@)), where t = 1 


3.3 Lottery Model 
In the lottery model, the temporal evolution of the morphogenetic field is trivial: 
fee =e 


The evolution of the geometric morphogenetic fields is then given by the initial conditions 


of the temporal evolution field: 


3.4 Advected Sandpile 


In the advected sandpile, the geometric morphogenetic fields are coupled with a dissipation 
field and given as: 


Oti t,k + (V fit + (V Nae 0 


Where f is the geometric morphogenetic potential, and h is a dissipation potential. This 
means that we can re-scale f and h by a constant factor Q to model the fact that they are 


both proportional to the rate of change of the fractal geometry of the evolving sandpile: 
Oii t,k + (V fit + (V -h)itk = O, 


Where k denotes the fractal dimension of the evolving sandpile. 


The dissipation field is written as w : R —> R, such that w(y) = 1 + +Q(y) with 
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\Q(y)| < 1. The dissipation is added to the advected sandpile to help keep the sandpile from 
collapsing because when it collapses, it will destroy the advective flux of morphogenetic 


potential through the sandpile. 


3.5 Dynamic Evolution 


It is now possible to write the dynamic evolution of the advected sandpile: 


It is worth noting that, if a # 0 (ie., if the advected sandpile does not collapse), then 
the above evolution equation will have a solution if and only if Q(x) < 1, where x = (x,y), 


namely, if and only if 


and Je > O(2,y) E RÉ = > |Q(x) - 1] < e. 


The first part of the above equation states that the sandpile has to expand away from any 
point in the domain (otherwise it will collapse) and then, under the dissipative term, the 
sandpile will shrink away from any other point and, finally, converge to any point on its own. 


These points are defined by 
{(x,y):0<a<A,y= fá- h} 
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This is the domain where the problem is expected to be most relevant and, as such, where 
it has to be solved by using the appropriate numerical scheme to solve the Navier-Stokes 
equation. Finally, the dissipative term is a conservative, linear, diffusion term that attempts 
to make the solution converge to the equilibrium state given by the initial conditions. It is 
worth mentioning that, in this case, the dissipation field is a random force applied on the 
advected sandpile and that the size of the sandpile is a random quantity with probability 


density given by 
p(r) = exp (- fw ar) 
0 


This has to be compared to the case of the advected sandpile: 


dp Fe exp (- I * 50) a) 


This shows the advantages of using a random force instead of the dissipative term. 


3.6 Numerical Implementation 


The numerical schemes developed in this paper are written using the advection and diffusion 
method, which allows the solution of the equation in time to be efficiently calculated, as 
it only requires knowledge of the flow velocity in the region around the advected sandpile 
at time tọ (so long as the solution of the advection equation is known). In this way, the 
numerical scheme is written using only vector and scalar fields. As an example, to write a 
linear scheme for the scalar advected sandpile, an explicit method for the evolution equation 


(1) is written in vector field representation as follows: 
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> 


q= So EA, Ho. 
For the linear scheme, we require the initial condition in vector field form, which can be 


written as: 


Ù = a Bo + 8 Ao. 


Here we have written the vector field By (and correspondingly Ap) as a linear combination 
of scalar fields, whose weights are related to the corresponding initial conditions. Using this 


linearization, the advection equations (1) can be written in the following matrix form: 


a 0 Blla 0 8 a 0 p 

0 a Bl 10 a 0 BL 10 a O B 
e olf olf o 

0 0 BY Os 0 8 00 6B 

0 0 0/0 0 0 0 0 0 


Using this matrix, the solution of the advection equation in time can be written as a 


single matrix-vector product with that matrix, giving: 
t = 5 — 7 Ho — a Bo + 8 Ay — a Bo + B Ay — y Ap — 8 Ao + Y Ao — B Ao — Y Ao + 8 Ao + 7 Ao 


We note that the solution of this equation is also the solution of the diffusion equation that 
arises when the sandpile is advected by a stream of fluid (with the same density distribution as 
the fluid and with no other specific properties; in other words, the sandpile moves according 
to a simple diffusion process). Similarly, the solution of the diffusion equation is also the 


solution of the equation (1) that would arise when the sandpile is advected by a stream of 
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fluid (with the same density distribution as the fluid and with no other specific properties; 
in other words, the sandpile moves according to a simple diffusion process). 
When the flow is laminar, we have that the matrix By vanishes, and we can simplify (1) 


by eliminating A: 


% = 75) — B Ao +7 Ao — 6 Ao + Ao + B Ao +7 Ao 


We find that, in this case, the solution of the diffusion equation becomes (for the sake of 
brevity): 

S(t) = A(t) — y Ap + 8 Ao 
where, by the third equation in (2), A(t) can be written as a linear combination of two inde- 


pendent functions (corresponding to the linear diffusion processes in the x-and y-directions): 
> T T 
A(t) = Ao + Ay exp(5¢) + Ao exp(5¢) 


The solution of (2) will then be the same solution of (1) that would arise when the sandpile 
is advected by a stream of fluid of density Ap (with no other specific properties) that moves 
according to the simple diffusion process that arises when the sandpile is advected by a 
stream of fluid of density Ag (with no other specific properties). However, in this case, the 
value of Ao is constant and equal to 1/m (where m is the mass of sand per unit area that 
is on the advected sandpile). Let us consider the case that the advected sandpile is initially 
flat. This situation is particularly interesting, because the solution of the problem (1) will 


show the behavior of the sandpile. In this case, we can use the definition of the stress tensor 


o= (At Ba, eA 
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and the definition of the stress tensor in the x-direction (that is, for the flow in the x-direction 
and for the x-direction alone): 


Org = 0 


We find that, in this case, the stress tensor is: 


= Ago! + Boor sale te (4t — Bo, | Ae — Aes 
= pan exP(an,) — Boor Yi exP( zm) — Boos 


Lx 


However, we can show that the stress tensor has the following solution (it is not unique, 


but a finite set of solutions): 


= Aost + Booz Boos l Aoor 1 
tt: j Ta 


— 2 2 (At — Bt)a, 


Since Ap has a finite norm (because it is bounded, hence positive), this is a solution of the 


problem. We also have that this solution satisfies the following relations 


OF x Z, Ao — 
Ta A Gila, 0) = oz(x,0), 


and 


Orz(£, +00) = 0. 


As a result, we find that the stress tensor o is finite everywhere (with a particular behavior 
near the sandpile boundary), but it can be decomposed in two components where the latter 
is a function of the shape of the sandpile and of the x-direction velocity. This is clearly a very 
intriguing problem to be analyzed. There are many interesting questions left unanswered, 


like for instance the behavior of the stress tensor near the x-direction boundary, the flow 
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induced by the sandpile, the formation of vortices, the influence of the advected sandpile on a 
gravity wave, etc. The solution that we have found for the stress tensor describes a situation 
in which the sandpile motion is forced to do horizontal water flows (with zero horizontal 
velocities) due to its own gravity and to a boundary force that acts as a sink of pressure (of 
the same order of magnitude as the sandpile density). Such a situation can be described by a 
parabolic equation (ø can be determined from the pressure drop and the horizontal velocity 
in the problem), and therefore it can be used to study the effect of the sandpile motion on 


the gravity wave field. 


4 Conclusions and Discussion 


Morphogenesis is an emergent phenomena that can only be described by the mathematics 
of nonlinear dynamics. This emergent phenomenon arises because of the discrete nature of 
spacetime and the continuous nature of physics. So far, all quantum gravity and most theo- 
retical approaches to modeling quantum gravity use a differential geometry, i.e., one where 
spacetime is an infinitely differentiable manifold. In these theories, the natural boundary 
conditions that should be used to find classical or semi-classical dynamics is that spacetime 
is asymptotically flat or anti-de-Sitter. However, these theories leave out the notion that 
the actual metric might take on some of these geometries, instead it is a constraint. This 
approach ignores the fact that the set of all of the possible geometries of the universe is an 
infinite set [14], so many of those geometries will give you a divergent (i.e., singular) metric. 


This divergence is due to the fact that the metric of spacetime is locally Minkowskian in 
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the smooth limit, but the overall structure of the spacetime geometry is locally fractal. This 
implies that all known theories of gravity, from general relativity, to string theory, to loop 
quantum gravity are local approximations of an underlying non-linear theory of gravity. This 
theory of gravity describes the dynamics of spacetime as a fractal system with a hierarchical 


set of symmetries, i.e., diffeomorphisms, i.e., coordinate transformations. 


A Fractal Dimension in the Anatomy of the Brain 


There have been attempts to apply fractal geometry to brain anatomy and function [4; 
2; 3; 6]. The fractal dimension of a geometric object is an exponent of the object’s box- 
counting fractal dimension, an approximation to its Hausdorff dimension. For example, 
the fractal dimension of a spheroid in 3-D is approximated by the exponent of the number 
of (nonoverlapping) spherical shells that fit inside its bounding box. So, an estimate of the 
(approximate) fractal dimension of a brain can be found by calculating the fractal dimension 
of the cortex as a whole using some appropriate mathematical framework for brain anatomy, 
and then finding the power law relating this "apparent" fractal dimension and the number 
of layers in the cortex. Fractal dimension then provides a dimensionless and continuous 
(though approximate) measure of brain structure. The idea may prove fruitful since it has 


recently been suggested that cortical dynamics might be fractal in nature [5]. 
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